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Two-dimensional non-abelian quantum field models provide a useful laboratory for analytic and 
numerical investigations of quantum theories with gauge symmetry. They can exhibit various fea- 
tures, such as charge confinement, which are known from D=4 theories like QCD. Several analytic 
predictions concerning the spectra of two-dimensional systems with adjoint matter were postu- 
lated and numerical results were obtained using Discrete Light Cone Quantization techniques, 
however none of them has been checked via Monte Carlo simulations. In this Letter we present 
two such models which are particularly interesting from the physical point of view and discuss 
first numerical results. 
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Grasping the full, nonperturbative dynamics of strong interactions is a difficult task. Although 
current cutting-edge lattice simulations of QCD provide the hadron spectrum obtained from first 
principles [j|], they do not unravel the hidden complexity of the theory. Hence, one would like to 
simplify QCD in order to disentangle some of the mechanisms at play and investigate them sepa- 
rately. From this point of view, two ways of approximating QCD appear as especially appealing. 
The first one is dimensional reduction, where one investigates lower dimensional theories instead 
of D = 4 ones. Although the former are much simpler, they may still possess some of intrinsic 
features of the original theory. The second way is the approximation which uses the large-N limit 
[0]. In this case, one replaces the SU (3) color symmetry group by SU (N) and studies the theory 
in the limit of N — > °°. An interplay of both these approximations enabled to construct several 
models having interesting properties. In the following we will describe two of them, which appear 
particularly attractive, and present some numerical results obtained for systems which are related, 
but technically simpler to simulate. 



1. Models of QCD 

1.1 Fundamental matter 

In 1962 Julian Schwinger solved analytically a two-dimensional version of quantum electro- 
dynamics [J3j] . He showed that massless electrons and positrons are confined into neutral pairs - 
massive mesons. In this model a free charge of any value is screened by pairs coming from the 
vacuum [Qj. Hence, the spectrum is composed exclusively of identical, massive mesons as is clear 
from the bosonised version of the model. Subsequently, generalizations for systems with Nf flavors 
of electrons were found [Q], in case of which one gets, with a nomenclature borrowed from QCD, 
a number of massless pions and a single massive eta meson. No exact solutions exist for models 
with massive electrons, however several approximate results were given [Q,|6|]. Non-abelian models 
were also considered []7|]. For the SU (3) gauge symmetry group, such system corresponds simply 
to dimensionally reduced QCD. 

In 1974 't Hooft proposed another theoretical tool to simplify QCD. He considered a two 
dimensional version of the theory in the limit of large number of colors [JsJ] and found that the 
complexity of the dynamics is drastically reduced when N is taken large. The resulting spectrum 
had several phenomenologically attractive features such as narrow, non-interacting mesons (see 
figure [[]). However, this approximation had also some drawbacks: quark loops were suppressed, 
hence rendering many body bound states invisible and problems with reconstruction of barionic 
spectrum appeared. 

From the numerical perspective a lot of work was devoted to the Schwinger model. On one 
hand, due to its technical simplicity it was used as a testbed for many numerical algorithms [Q, |10| ]. 
On the other hand, due to the relative complexity of its spectrum and of phenomena present in it, 
the Schwinger model gave a lot of understanding into the dynamics of confining gauge theories. 
Many complementary numerical results were obtained through different lattice calculations as well 



as with the discrete light cone quantization method (for the latter see for example Ref.[]ll[]). 
1.2 Adjoint matter 

It seems very likely that the interplay between a lower-dimensional theory and the large-N 
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Figure 1: Schematic plot showing the difference between bound states composed of fundamental (left) and 
adjoint (right) matter (dots) and glue fluxes (segments). Fundamental matter fields have a single color index; 
in the large N limit they can be bound in two body states only. Adjoint matter fields have two color indices; 
hence, many body bound states are possible. 



limit can lead to a model which can capture some essential features of QCD. A way to correct the 
drawbacks of the 't Hooft model is to change the representation of the matter fields. One can use 
quarks in the antisymmetric representation [12], which for N = 3 reduces to the fundamental one. 
Hence, one gets a theory which corresponds to QCD for N = 3 but possess a different large-N limit, 
i.e. quark loops are not negligible and the spectrum may contain many body bound states. Another 
possibility is to use adjoint matter fields. In this case, quark loops are of the same order as the 
gluonic ones in the expansion parameter 1 /N for N large. 

The two most common ways to introduce to the model matter fields in the adjoint representa- 
tion are either by dimensional reduction or by supersymmetry. In the former case one gets scalar 
fields as the components of the reduced gauge potential, whereas in the latter case one obtains 
fermions as the superpartners of the bosonic gauge degrees of freedom. Both possibilities were 
investigated in Refs.[|T3|]. 

The simplest model, which also attracted some attention, is composed of a single flavor of 
adjoint Majorana fermions [14]. Its spectrum contains many body bound states and an exponen- 
tially growing density of states - characteristic features of QCD spectrum. As far as the numerical 
results are concerned, only findings obtained by the discrete light front quantization technique are 
available [^] - no lattice simulations of this two dimensional model were performed up to now. 

Therefore, we set up a project of a systematic study by lattice methods of QCD2 models with 
adjoint Majorana fermions for different SU (N) gauge groups. Our aim is to crosscheck the discrete 
light cone quantization method results and verify theoretical predictions. We have chosen two 
particular models which were not studied on the lattice extensively and for which a connection to 
QCD can be argued. We briefly describe these two systems in the next section. 



2. Interesting models 
2.1 Kutasov model 

In Ref.[16] Kutasov considered a system consisting of a SU(N) gauge field and N 2 — 1 real 



Majorana fermions transforming in the adjoint representation of the gauge group, described by 
means of a traceless, hermitian N x N matrix. The spectrum is expected to consist of many body 
bound states (see figure [j]) with a density of states of a given mass raising exponentially with mass 
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[13]. The lagrangian is given by, 

N f 



2 = Tr(lF M 2 v + £ (ij^f-D^+mk^)). (2.1) 



jfc=i 

with g and m being two dimensionfull parameters of the dimension of mass. 

In two space-time dimensions there is no dynamical gauge degrees of freedom, hence the 
model contains only fermionic degrees of freedom. An analytic calculation in the light cone coor- 
dinate frame can reveal this. In the A'L = gauge the Lagrangian reads 

^ = Tr(^\(d-A + ) 2 + iyd + y + i\j/d-\jf-2im\jf\jf+A + J + ^, (2.2) 

where an d %j denote respectively the right and left moving fermions. After integrating the non 
dynamical degrees of freedom one gets 

JSf = Tr (iyfd+ iff + g 2 J + - im 2 y/-^- (2.3) 

which is exclusively composed of right-moving fermions. 

On first sight the system eq.(23) cannot be supersymmetric, since the numbers of bosonic and 
fermionic degrees of freedom. Nevertheless, it turns out that for a particular value of the quotient 
g/m the model possess a dynamical supersymmetry [16]. Namely, there exist two supersymmetry 
generators given by 

1 f°° 

Q + = —= \ dp dq\i/ij(-p)Wjk(q)¥kiip ~ q), 

Q-=g[ ^-^-Yij{-p)Yjk(q)Yki(p ~ q), (2.4) 

J -co p 

satisfying for m 2 = g 2 N 

[Q+,p-]=0, [Q-,p-]=0, 

{Q + ) 2 =P + , (Q-) 2 = -P-, (2-5) 
P and P + being the light-cone Hamiltonian and total momentum operators respectively. Moreover 

{Q + ,Q~}=0. (2.6) 
Finally, one should note a result published in Ref. Jl7[] pointing to a partial equivalence between 



the model eq.(2.1) and a model with fundamental matter. Namely, the Authors proved that the 
massive sector of a system with SU (N) color symmetry and a single adjoint fermion is equivalent to 
the massive sector of a system with SU (N) color symmetry and N flavors of fundamental fermions. 

2.2 = (1,1) model 

The second model of interest is a two-dimensional theory obtained by dimensional reduction 
of the D = 4, ,jV = 1 supersymmetric Yang-Mills theory down to D = 2. It was considered as a 
toy model of QCD by Dorigoni, Wosiek and Veneziano [18], where it was analyzed numerically 
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Figure 2: Pion mass and fermion condensate for the Schwinger model, comparison with Ref . pOj] (parame- 
ters: 8x20 lattice, 10 4 measurements). 



by the light-cone quantization method. This system was also studied on the lattice in Ref.fJI^] in 
the context of supersymmetry on the lattice. 

The D = 4 theory reduced to D = 2 yields a lagrangian of the following form, 

Se = Tr ( - -F^F^ + /AY D^A + i% 

+ 2 + \ D v & Dli & + mteration terms) . (2.7) 

Hence, it may be viewed as a generalization of the Kutasov model by inclusion of a second fermion 
and two scalar fields. 

Using the numerical light-cone approach and the Coulomb approximation the model was 
shown to be a non-trivial generalization of the 't Hooft's model with an arbitrary number of par- 
tons. The spectrum turned out to be composed of color-less bound states of many partons whose 
mass grows linearly with distance between the partons, namely M ~ a£ | Ax| [|T8|]. Hence, a string 
picture of flux tubes connecting each parton emerged. 



3. First results and plots 

We have started our project by setting up simulations with two Dirac adjoint fermions. We 
have implemented dynamical fundamental and adjoint fermions using the Hybrid Monte Carlo al- 
gorithm. A first set of tests consisted in reproducing the known results for the Schwinger model. On 
figure |2| we compare our results with the numbers published in Ref.[2C]. A very good agreement is 



seen. Then, we have performed simulations for the SU (2) gauge group. Figure || show the compar- 
ison of the pion mass for the U(l), SU (2) fundamental and SU (2) adjoint models. The equivalence 
of Ref.JT/]] and conclusions of Ref. [|5[] suggests that the massive states of the SU (2) adjoint model 
should be \/2 times heavier than the massive states of the SU (2) fundamental model. Data shown 
on figure || indeed hint that the adjoint pions are heavier than the fundamental ones. However for 
the exact determination of the coefficient studies in the v\ sector are needed and extrapolations to 
the continuum and to the chiral limit must be performed. 
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Figure 3: Pion masses for f/(l) and SU(2) fundamental and adjoint models (parameters: 8x20 lattice, 10 4 
measurements). 

4. Conclusions 

Summarizing, we argued that two-dimensional Yang-Mills theories enable to grasp particular 
features of quantum chromodynamics. By appropriately choosing the representation of fermionic 
matter fields it is possible to obtain models with an interesting large-A^ limit. Through systematic 
studies one may find a system which resembles QCD particularly well, but thanks to its simpli- 
fied structure may be analytically tractable and lead to a better understanding of the mechanisms 
at play. Lattice approach is particularly well suited for such investigations. We recalled two two- 
dimensional models, a supersymmetric Yang-Mills theory obtained by dimensional reduction, and 
a non-supersymmetric Yang-Mills theory with adjoint matter fields which turns out to be super- 
symmetric for certain value of parameters, as particularly interesting to study numerically. Finally, 
we discussed first numerical results. Further studies are on the way. 
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